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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2213 

AERODYNAMIC  COEFFICIENTS  FOR  AN  OSCILLATING  AIRFOIL  WITH 
HINGED  FLAP,  WITH  TABLES  FOR  A  MACH  NUMBER  OF  0.7 
By  M.  J.  Turner  and  S.  Rabinowitz 

SUMMARY 


Dietze*s  method  for  the  solution  of  Possio*s  integral  equation  has 
been  used  to  detemine  the  chordwise  distribution  of  lift  on  an  oscil¬ 
lating  airfoil  with  simple  hinged  flap  in  two-dimensional  compressible 
flow  (subsonic).  The  results  of  these  calculations  have  been  used  to 
prepare  tables  of  aerodynamic  coefficients  for  lift,  pitching  moment 
(referred  to  quarter-chord  point),  and  flap  hinge  moment  for  a  Mach 
number  of  0.7;  the  motions  considered  are  vertical  translation,  airfoil 
rotation  about  the  quarter-chord  point,  and  rotation  of  the  flap  about 
its  hinge  line. 

Aerodynamic  coefficients  are  tabulated  for  k  values  of  (ratio 

of  flap  chord  to  total  chord)  and  for  12  values  of  reduced  frequency 
covering  the  range  from  0  to  0.7.  Results  are  given  for  one  value  of 
Mach  number,  M  =  0.7.  Data  of  this  kind  have  already  been  presented 
by  Dietze  for  one  value  of  the  ratio  of  flap  chord  to  total  chord, 

=  0.15;  these  results  have  been  checked  independently,  and  calcula¬ 
tions  have  been  carried  out  for  three  additional  values,  =  0.24, 

0.33>  and  0.42.  Certain  auxiliary  parameters,  which  will  be  needed  in 
any  further  calculations  of  this  type,  are  presented  for  future  reference. 


INTRODUCTION 


The . fundamental  integral" equation  for  the  pressure  distribution  on 
an  oscillating  thin  airfoil  moving  at  subsonic  speed  has  been  derived 
by  Possio  in  reference  1.  Collocation  procedures  have  been  used  by 
Possio,  Frazer  and  Skan,  and  others  to  obtain  lift  and  moment  on  an 
oscillating  flat  plate.  An  important  contribution  has  been  made  by 
Dietze  (see  references  2  and  3)^  who  has  developed  an  iterative  procedure 
for  numerical  solution  of  Possio *s  integral  equation.  This  procedure  is 
particularly  well  adapted  to  the  calculation  of  aerodynamic  loading  on 
an  oscillating  airfoil  with  hinged  flap. 
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It  has  been  pointed  out  correctly  by  Karp  and  Weil  (reference  h, 
p.  ll)  that  Dietze's  procedure  does  not  properly  account  for  the  loga¬ 
rithmic  singularity  in  pressure  distribution  at  the  flap  hinge.  However, 
for  applications  in  flutter  analysis  the  principal  objective  is  to  deter¬ 
mine  resultant  lift  and  moments  rather  than  pressure  distribution.  Mathe¬ 
matically  exact  solutions  in  closed  foim  are  available  for  the  stationary 
thin  airfoil  with  deflected  flap,  and  it  is  found  that  lift  and  moments 
obtained  by  Dietze's  iterative  procedure  are  in  excellent  agreement  with 
theoretically  exact  values.  There  is  good  reason  to  expect  that  equally 
satisfactory  results  will  be  obtained  for  the  airfoil  with  oscillating 
flap.  The  existence  of  the  singularity  in  pressure  distribution  is  a 
consequence  of  the  shaip  corner  in  the  idealized,  broken-line  profile, 
and  of  the  associated  discontinuity  in  downwash  velocity  at  the  flap 
hinge.  The  introduction  of  a  finite  cosine  series  for  the  downwash  is 
in  effect  equivalent  to  a  slight  modification  of  the  idealized  profile 
by  rounding  off  the  comer  at  the  flap  hinge. 

This  work  has  been  performed  at  Chance  Vought  Aircraft,  under  the 
sponsorship  of  the  Bureau  of  Aeronautics,  Navy  Department,  in  order  to 
provide  data  for  the  calculation  of  con5)ressibility  effects  in  control- 
surface  flutter  problems.  It  has  been  made  available  to  the  National 
Advisory  Committee  for  Aeronautics  for  publication  because  of  its  general 
interest. 


SYMBOLS 


ratio  of  flap  chord  length  to  total  chord  length 


Ra 


airfoil  region  in  x,z-plane 


5y(x,t)  vertical  displacement  of  point  on  idealized  profile, 

positive  upward 

6(1)  nondimens ional  representation  of  instantaneous  chordwise 

distribution  of  vertical  displacement 

I  total  chord  length 

03  circular  frequency 

t  time 

X  chordwise  coordinate 
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y 

z 

I 

Vy(x,t) 

g(|) 


V 

7(1) 

^inc 

K(s,M) 

u,v 

s 

M 

a  =  1  -  \fi 

p 

a(  x,t) 

T(av) 

AK(s,M) 

AK3^(s,M) 

AK2(s,M) 

kgn 


vertical  coordinate 
spanwise  coordinate 

dimensionless  chordvise  coordinate  (2x/2) 

vertical  component  of  fluid  velocity  adjacent  to  airfoil 

function  representing  instantaneous  distribution  of 
vertical  fluid  velocity  adjacent  to  airfoil 

(Vy(x,t)  =  8(1)6^“*) 
reduced  frequency  (cnZ/2V) 
velocity  of  flight 

function  representing  distribution  of  dipole  lines 
distribution  of  dipole  lines  for  incompressible  flow 


kernel  of  Possio^s  integral  equation 
variables  of  integration 
auxiliary  variable  ^0)^(5  -  ^o)) 


Mach  number 


air  density 
lift  per  unit  area 

function  defined  by  Kussner  and  Schwarz 

kernel  difference  ^K(s,M)  -  K(s,0)) 

singular  part  of  kemel  difference 

nonsingular  part  of  kemel  difference 

constant  occurring  in  formula  for  AK-^ 

coefficient  in  polynominal  representation  for  AK2 


k 
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Pik>  ^nv 


AMp 

AMr 


Cgh^  ^gh 


coefficients  in  recursion  formulas  for  solution  of 
Possio’s  integral  equation 

lift  force  on  airfoil  strip  of  width  Az 

pitching  moment  on  airfoil  strip  of  width  Az 

hinge  moment  on  flap  for  strip  of  width  Az 

downward  displacement  of  quarter-chord  point  divided  by 
semi  chord 

rotation  of  airfoil,  positive  in  stalling  direction 
rotation  of  flap,  positive  in  stalling  direction 
aerodynamic  coefficients,  where  g,  h  =  S,  D,  R 


BASIC  THEORY 


A  very  complete  digest  of  the  literature  on  aerodynamic  theory  of 
oscillating  airfoils  has  been  presented  by  Karp,  Shu,  and  Weil  in  refer¬ 
ence  5.  Consequently,  the  basic  theory  is  outlined  briefly  merely  to 
exhibit  the  essential  features  of  the  computational  scheme  and  to  point 
out  certain  errors  which  have  been  discovered  in  Dietze^s  formulas. 

The  usual  assuir5)tions  of  thin-airfoil  theory  are  adopted,  leading 
to  Possio*s  integral  equation  for  the  chordwise  lift  distribution  on 
the  oscillating  airfoil.  Rectangular  coordinates  are  employed,  with 
the  x-axis  alined  in  the  direction  of  the  undisturbed  flow.  The  airfoil 
is  replaced  by  a  deformable  sheet  of  zero  thickness  which,  in  its  undis¬ 
turbed  position,  occupies  the  region  Rg^ 

-  2/2  =  X  =  2/2,  y  =  0 


of  the  x,z-plane. 

The  lifting  surface  executes  sinusoidal  oscillations  in  which  each 
point  moves  along  a  line  parallel  to  the  vertical  y-axis.  Displacements 
are  independent  of  the  spanwise  coordinate  z,  and  may  he  represented  in 
the  form 


iCDt 


^y  “  "  2 
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In  accordance  with  the  usual  convention,  it  is  the  real  part  of  equa¬ 
tion  (l)  which  has  physical  significance.. 

The  y-component  of  velocity  of  a  fluid  particle  adjacent  to  the 
lifting  surface  is  related  to  5y  by  the  equation 


Vir  = 


Sby 

sr 


V 

OX 


5&y 

'ST 


— 

i 


5&y 


(2) 


where  i  =  2x/Z  and 
(2)  it  follows  that 


1=1  =  1  inside  Rg^.  From  equations  (l)  and 


Vy  =  gU)e^^ 


(3) 


where 


g(|)  =  vlio^y 


w 


where 


cOp  =  CJDZ/2V 


In  Rietze^s  derivation  of  the  "basic  integral  equation  the  region  Rg^ 
is  covered  with  dipole  lines  of  density  V7(x)ei“’*^  per  unit  length  in 
the  chordwise  direction.  By  equating  the  vertical  velocity  induced  by 
the  dipole  covering  to  that  given  by  equation  (3)  the  following  integral 
equation  is  obtained  for  the  determination  of  y: 

g(|)  7(Io)K(s,M)  d^o  (5) 

subject  to  the  condition  that  7(1)  shall  be  finite;  the  kernel  of  the 
integral  equation  is  given  by 


K{s,M)  =- 


isXJH 


k\jl  -  m2 


Hq^^HisIX,)  -  xM  |s|x) 

I  S  I 


.  X 

id  -  M2)e‘''®M 


log 


M 


jt  \/i  -  m2  1  -  \Ji  - 


eiu^Q(2)(|u|j4^  du 

Jo 


(6) 
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with 

s  =  -  ^o) 

\  =  M/(l  -  m2) 

The  lift  a  per  unit  area  (positive  upward)  is  given  hy 

a(x,t)  =  pV7(x)e^^'  (?) 

NUMERICAL  SOLUTION  OF  POSSIO^S  INTEGRAL  EQUATION 


Dietze's  approximate  solution  of  equation  (5)  (Possio*s  integral 
equation)  is  of  the  form 


^  ^  ^inc  +  7i  +  72  +  •  •  •  +  7n 


(8) 


where  7^^^  and  7^  { v  =  1,  2,  .  .  . ,  n)  are  solutions  of  the  integral 

equations 


g(l)  =  (Or  ^  7lnc(^o)K(s,0)  d|o 

-1 


(9) 


gyCO  =  <^^0^  V  =  2,  .  .  n  (10) 

-1 


and  where 


Ol 


gl(|)  =03^'?)  7inc(^o)[K(s,0)  -  K{s,M)]  d|o 

J-1 


(11) 


nl 

gyd)  7v-i(Io)[k(s,0)  -  K(s,M)]  d|o,  V  =  2,  3,  •  •  n 

J-1 


(12) 
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K(s,0)  =  lim  K(s,M) 

M~>0 

The  required  solutions  of  equations  (9)  and  (lO)  are  obtained  by  the 
methods  of  reference  6.  Convergence  of  the  process  has  been  proved  only 
for  the  stationary  case.  However,  computational  experience  furnishes 
convincing  evidence  that  the  process  does  converge  in  the  more  general 
case. 

It  will  be  observed  that  Dietze's  process  requires  essentially  the 
solution  of  a  succession  of  integral  equations  with  kernel  K(s,0)  from 
the  incon^ressible  problem.  The  functions  g^(|)  are  obtained  by 
direct  integration  in  accordance  with  equations  (ll)  and  (12). 

In  case  can  be  represented  by  a  cosine  series  of  form 

gv(l)  =  V^O  cos  nj^Y  0  =  (Zf  =  jT  (I3) 

\  n=l  / 


with 


I  =  -cos  0 


then  it  is  known  from  the  work  of  Kiissner  and  Schwarz  (reference  6)  that 

=  -2v|ao  cot  ^  +  2  ®n  (l^) 


where 


-  (  5 — \^0  -  -^1 )  + 


^  )  -  An,  n  =  1 


and  T(cDr)  is  the  function  of  reduced  frequency  defined  in  reference  6. 
The  T-function  is  related  to  Theodorsen's  C-function  by  the  equa¬ 
tion  T  =  2C  -  1. 


In  order  to  facilitate  the  evaluation  of  the  integral  occurring  in 
equation  (l2)  the  kernel  difference 


AK(s,M)  =  K(s,M)  -  K(s,0) 
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is  expressed  in  "the  form 

AK(s,M)  =  Z^]^(s,M)  +AK2(s,M) 

where 


AKi(s,M)  =  ^  +  kii  +  ki2  loge  |s|  +  s^ki3  +  ki4  loge  |s|)  (l?) 


loge 7  =  0.57T22  (Euler's  constant) 
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The  nonsingular  part  AK2(s,M)  is  replaced  by  a  polynomial  of  ninth 
degree-*- 


AK2(s,M)  «  -  kpns** 


whose  coefficients  are  determined  separately  for  each  Mach  number  by 
fitting  the  polynomial  to  tabulated  values  in  the  interval  |s|  =  1.8. 

Equation  (l2)  may  be  written  in  the  form 


gy(l)  =  7y_j  AK(s,M)  d^o 


If  it  be  assumed  that 


7y_i  =  -2v/po  cot  I  +  2  XI  Pn  sin 


then  it  follows  that  (upon  carrying  out  the  required  integrals 
defining  g^,  in  accordance  with  equations  (1T),  (18),  and  (l9)  and 
applying  equations  (15)  to  solve  equation  (lO)), 


=  -2V^c 


where 


9n  -  9in  +  92n 


^ietze's  definition  of  the  approximating  polynomial  (table  3  of 
reference  2)  differs  in  sign  from  that  given  here.  However,  it  has  been 
found  by  carefully  checking  the  derivations  that  the  minus  sign  is 
required  in  equation  (l8)  in  order  to  justify  the  recursion  formulas 
for  calculating  7v  from  It  is  believed  that  this  is  merely  an 

error  of  presentation,  since  Dietze*s  numerical  results  are  found  to  be 
substantially  correct. 
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Presmably  this  is  merely  an  error  in  presentation,  since,  as  already 
noted,  Dietze*s  numerical  results  have  been  verified  by  independent 
calculations.  Numerical  values  of  the  parameters  a^,  for  M  =  0.7 

and  reduced  frequencies  ranging  from  0  to  O.7  are  given  in  table  I. 


EVALUATION  OF  KEBHEL  AND  KERNEL  DIFFERENCE 


By  making  use  of  the  relations  among  Bessel,  Hankel,  and  Neumann 
functions  and  by  separating  the  kernel  into  real  and  imaginary  parts  an 
ejq)ression  of  the  following  form  is  obtained: 

K(s,M)  =  K'(s,M)  +  lK"(s,M)  =  cos  (sXM)  +  B(s,M)  sin  (sXM)  ^ 


h\Jl  -  m2 


i  B(s,M)  cos  (sAJ^)  -  A(s,M)  sin  (sAM) 

k\Jl  -  m2 


(25) 


where 


A  -  Jq(  I  s  W)  -  M  N2^(  I  s  I  X)  -  sjl  -  log^ 


M 


1  -  \}l  -  M‘ 


sin 


H)- 


(1  -  m2)  (“m).L  [cos  u  Jo(|u|M)  +  sin  u  No(|u|M)]  du 


4 


(1  -  m2)  cos  ^s  J  [cos  u  Ho(|u|M)  -  sin  u  JqC |u|M)]  du 
B  ^  Nod  six)  +  M-rfr  Ji(lslX)  +  §\/l  -  m2  loge  - ^  cos  fs 


(26) 


M 


(1  -  M^)  cos  ^s  [cos  u  Jo(|ulM)  +  sin  u  No(|u|M)J  du  - 

0 


X 

IS— 

M 


(1  -  m2)  sin  fs  Qcos  u  No(  |u(M)  -  sin  u  Jo(  lu|M^  du 
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By  introducing  a  new  variable  of  integration 


M 

"  X  ^ 


equations  (26)  are  transformed  into 


A  =  Jo(|s|X)  -  M  jIj  NidslX)  -  |\/l  -  m2  loge 


1  -  Vl  -  m2 


G  m)  - 


(s  j  jcos  1)  Jo(|v|X)  +  sin  |)no(  |v|X)J  dv  - 


COS  IS 


m) 

t  0  ^ 


|v|X.)  -  sin  (v 


i)jo(lvU)l 


lIo(  1»|X)  +  M  ^  Jid.l).)  +  fVl  -  lose  ^  ^  OO'  ('  h)* 

iS 

cos  ^s  j^cos  (v  1)  Jo(lv|X)  +  sin  (v  |)no(  |vIX)J  dv  - 

sin  (s  J  ^os  (v  No(  |v|X)  -  sin  (v  |)  JqC  |v|  X)j  dv  (2{ 
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Numerical  values  of  the  following  Integrals  are  required: 


Il(s)  = 


COS 


l2(s)  =  j  sin  Jo(|v|X)  dv 


13(3)  = 


I4(s) 


TS 

do 

ns 


cos 


sin 


(v  No{  |v|X)  dv 
^v^)no(|v|X)  dv 


(29) 


The  evaluation  of  and  can  be  obtained  by  numerical  Integration 

In  a  stral^tforward  manner.  However^  since  No(x)  has  a  logarithmic 
singularity  at  x  =  0,  It  Is  necessary  to  express  and  In 

different  form.  Upon  making  the  substitution  (reference  PP.  130,  132) 

§  NqCx)  =  Jo(x)  loge  ^  -  Bo(x) 


where 


Bn(x)  =- 


(If  *  d(i)'  - 


,  1  1 
(3.')2 


(Ip  ■  • 


Ik 
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and  in'tegra’ting  by  parts,  the  following  equations  ane  obtained: 


I  13(5)  =  Ii(s)  loge  ^ 
J  14(3)  =  l2(s)‘  loge  ^ 


It  follows  from  equations  (29)  that  Ii  and  I3  are  odd  functions, 
while  I2  and  are  even.  The  integrals  occiirring  in  equations  (30) 

and  ( 31)  can  be  evaluated  numerically  without  difficulty;  it  should  be 
noted  that 


lim 

v-->0 


Il(v) 

V 


=  1 


lim 

v->0 


0 


In  computing  numerical  values  of  the  kernel  Dietze  has  used  the 
tables  of  Bessel  and  Neumann  functions  given  in  reference  7,  which  do 
not  permit  a  satisfactory  determination  of  Bq  since  Nq  is  tab\ilated 
to  only  four  (in  some  cases  three)  places.  In  recalculating  the  kernel 
the  seven-place  tables  given  in  reference  8  have  been  used. 

For  evaluation  of  the  integrals  Iq  the  form\ilas  given  in  refer¬ 
ence  9,  page  227,  have  been  extended  to  include  fifth  differences.  These 
formulas  have  been  used  with  an  interval  A(sX)  =  0.05,  and  the  results 
have  been  checked  up  to  sX  =  O.3O  by  using  an  interval  A(sX)  =  0.02. 
Recalculated  values  of  K(s,0.7),  K(s,0),  AK(s,0.7),  AK2^(s,0.7),  and 
AK2(s,0.7)  are  given  in  table  II.  These  values  are  found  to  be  in 
close  agreement  with  those  given  by  Dietze;  where  differences  exist,  the 
new  values  are  believed  to  be  more  accinrate.  In  making  comparisons  it 
should  be  noted  that  Dietze  has  tabulated  -K(s,M). 

The  coefficients  k2n  are  obtained  by  fitting  a  ninth-degree 
polynomial  (see  equation  (I8))  to  the  tabulated  values  of  AK2  in  the 


MCA  TN  2213 


15 


interval  |sl  =  1.8  by  the  method  of  least  squares*  The  values  obtained 
in  this  way  for  M  =  O.7  are 

k22  =  -0.0lf6728  -  O.Olt.5974i 

k23  =  0.023019  -  0.023^911 

k2l^  =  0.009818  +  0.0528551 

k25  =  -0.020228  +  0. 003^1 

k26  =  -0.001040  -  0. 0215101 

k2Y  =  0.007272  -  0.0002831 
k23  =  0.000053  +  0.0031171 
k29  =  -0.000997  +  0.0000121 

Since  |s|  =  -  Iq)  ||  -  |  =  2  the  approximate  representa¬ 

tion  for  AK2  is  valid  for  reduced  frequencies  in  the  range  0  =  0)1.  =  0.9. 
Although  the  values  of  the  coefficients  k2n  given  here  differ  from  those 
given  by  Dietze  (see  translation  of  reference  2,  p.  26),  there  is  reason¬ 
able  agreement  of  coefficients  for  lower  powers  of  s.  Also  the  algebraic 
signs  agree;  this  gives  further  indication  that  Dietze  must  have  intro¬ 
duced  a  minus  sign  in  his  definition  of  the  approximate  representation 
for  AK2 . 


NOTATION  FOR  AERODYMMIC  LIFT  AND  MOIMTS 


In  presenting  his  nmerlcal  results  Dietze  has  Introduced  repre¬ 
sentations  of  lift  force  and  moments  involving  only  real  quantities. 
Complex  notation  is  used  in  the  present  report  in  order  to  conform  more 
nearly  to  ctirrent  American  practice;  however  the  essential  features  of 
Dietze 's  notation  (translation  of  reference  3)  are  retained.  Lift  and 
moments  on  an  airfoil  strip  of  width  Az  for  an  airfoil  with  simple 
hinged  flap  are  as  follows: 


l6 
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APg  =  Az|^t0^2cgg  _  kgg^qg  +  (‘Oj.^cgj)  -  kg^^qj,  +  (o^r^cgj,  -  kgR)q^ 

AMj)  =  -  kjjg^qg  +  (tUp^Cjio  -  kjjp^qp  +  -  kjjp^qpJ 

AMg  =  itpv2  Az|(“/cgg  -  kpg^qg  +  (uj/cpp  -  kgp^qp  +  -  kp^^q^J 

where 


APs 

AMj) 


AMr 


lift  force,  positive  upward 

stalling  moment  on  airfoil  plus  flap,  referred  to  quarter- 
chord  point 

hinge  moment  on  flap,  positive  in  same  sense  as  AMp 


©' 


downward  displacement  of  quarter-chord  point 

rotation  of  airfoil  in  stalling  direction 
rotation  of  flap  in  stalling  direction 


kgh  =  l^gh*  + 


qh  =  (qh’  +  iqh")®^* 


g,  h  =  S,  D,  R 


The  quantities  Cgj^  may  he  expressed  in  terms  of  the  functions  $i^, 
<S>rjf  and  $2^2  defined  in  reference  6  as  follows: 


S 

D 

R 

s 

1 

1/2 

$4/2m 

V 

1/2 

3/8 

R 

$l(,/2jt 

0  y/^TC 

.  . 

H32) 
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DISCUSSION  OF  NUMEEICAL  RESULTS 

The  coefficients  ^33,  ^DS^  ^SD^  ^DD^  which  do  not  depend  on 

the  ratio  of  flap  chord  to  total  chord  are  presented  in  table  III 

and  in  figures  1  to  8  for  M  =  O.7  and  a  range  of  from  0  to  0.7* 

The  hinge-moment  coefficients  and  kj^  associated  with  air¬ 

foil  flapping  and  rotational  motions  are  presented  in  table  TV  and 
figures  9  to  12  for  M  =  0.7^  reduced  frequencies  from  odj.  =  0  to  0.7^ 
and  ratios  tr  =  0.15,  0.24,  0-33,  an^  0.42.  Coefficients  in  the 

series  representation  for  7 


are  presented  in  table  V  for  both  flapping  and  rotational  motions.  It 
is  noted  that  the  coefficients  kR3  and  kj^  may  be  computed  for  any 

ratio  of  flap  chord  to  total  chord  without  further  iterations  by 
inserting  the  coefficients  for  the  appropriate  type  of  airfoil 

motion  into  the  formula 

kph  =  (cRhO)/  h  =  S,  D 

The  quantities  are  expressed  as  functions  of  Tjj  as  follows: 

cos  6  =  2tr  -  1,  0  9  ^  IT 


Qo  = 

(jT  - 

e)i 

-1+2  cos 

e) 

+  2 

sin  6  - 

-  sin  20 

2 

Ql  = 

-  e) 

COS  6  + 

3 

2 

sin  0 

1 

sin 

36 

Q2  = 

-(jt 

-  e) 

2  . 
--sin 

29 

+ 

1 

12 

sin 

n  - 

sin 

(n 

-  2)0 

2 

sin 

n0 

-1- 

sin  (n  + 

2)0 

Hn  - 

(n  - 

l)(l 

a  -  2) 

(n 

.  - 

1)(: 

n  + 

1) 

(n  +  l)(n 

+  2)' 

n  >  2 
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The  aerodynamic  coefficients  associated  with  rotational  motion  of 
the  flap  ^31^,  ^ER  presented  in  table  VI  and  figures  13 

to  l8  for  M  =  0.7,  reduced  frequencies  ranging  from  (Dj.  =  0  to  0.7^  and 
ratios  of  flap  chord  to  total  chord  -  0.15,  0.24,  0.33,  and  0.42. 

From  five  to  eight  iterations  have  been  employed  in  the  calculation 
of  7.  In  computing  the  coefficient  kj^  the  accuracy  depends  on  the 

number  of  coefficients  used  as  a  starting  basis  for  The  number  of 

coefficients  is  reduced  by  3  in  each  successive  iteration;  in  the  calcu¬ 
lations  described  in  this  report  22  coefficients  have  been  used  as  a 
starting  basis.  All  coefficients  available  at  a  given  stage  of  the 
iteration  have  been  used  to  calculate  the  contribution  of  to  7 

and  the  contribution  of  7±nc  been  obtained  in  closed  form  from  the 

results  of  reference  6. 

In  recalculating  the  coefficients  which  are  independent  of  and 

of  the  remaining  coefficients  for  =  0.15,  values  have  been  obtained 

which  differ  in  general  by  less  than  1  percent  from  those  given  by 
Dietze.  There  are  a  few  isolated  exceptions,  however.  An  error  of 
7  percent  has  been  foimd  in  the  imaginary  part  of  kj)3  for  cDp  =  0.10. 

Also  there  are  errors  in  Dietze  *s  values  of  the  imaginary  part  of  kj^ 

at  oij*  =  0.02  for  all  Mach  numbers  tabulated,  including  M  =  0.  The 
entry  for  M  =  0  should  read  lO^kj^  ^  ^  =  I.I8  instead  of  I.03. 

Apparently  the  same  error  has  been  carried  through  for  all  Mach  numbers. 

A  similar  error  occurs  in  the  imaginary  part  of  kj^  at  oip  =  O.6O;  the 

entry  for  M  =  0  should  be  lO^j^^*  =  131-8  instead  of  132.8,  and 

this  error  has  been  carried  through  for  other  values  of  Mach  number  as 
well . 


Chance  Vought  Aircraft 

Division  of  United  Aircraft  Coipio  rat  ion 
Dallas,  Tex.,  July  19,  1949 
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TABLE  I.-  VALUES  OF  AMD  FOR  M  =  0.? 

j-ai  =  ai'  +  iai'-j  %k  =  3ik’  +  iPik"] 


o-l  “2  “3 


CDj. 

n  *  * 

02* 

02  *  * 

«'3" 

4 - 

0.02 

0 

-0.04801 

-0.00094 

0 

0 

i  0.00001 

.04 

0 

-.09601 

-.00376 

0 

0 

.00006 

.06 

0 

-.14402 

-.00847 

0 

0 

.00021 

.08 

0 

-.19202 

-.01506 

0 

0 

.00049 

.10 

0 

-.24003 

-.02353 

0 

0 

.00095 

.20 

0 

-.48006 

-.09418 

0 

0 

.00762 

.30 

0 

-.72009 

-.21179 

0 

0 

.02573 

.40 

0 

-.96012 

-.37652 

0 

0 

.06099 

.50 

0 

-1.20015 

-.58831 

0 

0 

.11912 

.60 

0 

-1.44018 

-.84716 

0 

0 

.20584 

.70 

0 

-1.68021 

-1.15308 

0 

0 

.32686 

ox. 

Poo 

Po 

1 

Poo’ 

Poo" 

Poi' 

Poi" 

0.02  -1.82233  o.4o4i5  -0.00347  0.00169 

,o4  -1.63609  .64055  -.01064  .00688 

.06  -1.45768  .80639  -.01952  .01519 

.08  -1.29103  .92767  -.02917  .02624 

.10  -1.13670  1.01832  -.03904  .03961 

.20  -.52194  1.23247  -.08348  .13074 

.30  -.0779^<-  1.27472  -.11434  .24736 

.40  .27370  1.24893  -.13107  .38037 

.50  .57113  1.18299  -.13^19  .52669 

.60  .83322  1.08608  -.12388  .68538 

.70  1.06903  .96176  -.09990  .85624 
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TABLE  I.-  VALUES  OF  ai  ATO  Pik 
FOR  M  =  0.7  -  Continued 


^02 

^10 

P02’ 

P02" 

Pio' 

Pio" 

0.02 

-0.00001 

-0.00001 

-o.oo4i6 

0.00148 

.Oi^ 

-.00004 

-.00009 

'  -.01401 

.00591 

.06 

-.00018 

-.00024 

-.02809 

.01331 

.08 

-.00042 

-.00048 

-.04561 

.02366 

.10 

-.00079 

-.00081 

-.06601 

.03696 

.20 

-.00519 

-.00360 

-.19879 

.14786 

.30 

-.01438 

-.00769 

- . 36140 

.33268 

.40 

-.02896 

-.01257 

-.53416 

.591^3 

.50 

-.04790 

-.01796 

-.70335 

.92411 

.60 

-.07225 

-.02363 

-.85837 

1.33071 

.70 

-.10161 

-.02952 

-.99061 

1.81125 

o^r 

311 

3i2 

Pll' 

Pll" 

P12' 

P12” 

0.02 

-0.00001 

-0.04802 

0.00047 

0 

.04 

-.00005 

-.09611 

.00188 

0 

.06 

-.00016 

-.14430 

.00424 

0 

.08 

-.00038 

-.19262 

.00753 

0 

.10 

-.00075 

-.24108 

.01177 

0 

.20 

-.00599 

-.48583 

.04706 

0 

.30 

-.02020 

-.73434 

.10590 

0 

.40 

-.04790 

-.98513 

.18826 

0 

.50 

-.09356 

-1.23571 

.29416 

0 

,60 

-.16166 

-1.48287 

. 42358 

0 

.70 

-.25672 

-1.72280 

.57654 

0 

22 


mCA  TN  2213 


TABLE  I.-  VALUES  OF  ADD 

FOR  M  =  0.7  -  Concluded 


Pl3 

P20 

r 

Pis’ 

Pl3” 

0 

CVJ 

oa 

P20" 

0.02 

0 

-0.000004 

0.000003 

0.00001 

.04 

0 

-.00003 

.00002 

.00005 

.06 

0 

-.00010 

.00008 

.00014 

.08 

0 

-.00024 

.00019 

.00030 

.10 

0 

-.00048 

.00037 

.00053 

.20 

0 

-.00381 

.00299 

.00288 

.30 

0 

-.01286 

.01010 

.00713 

.1+0 

0 

-.03049 

.02395 

.01251 

.50 

0 

-.05956 

.04678 

.01778 

.60 

0 

-.10292 

.08083 

.02134 

.70 

0 

-.16343 

.12836 

.02130 

P21 

^22 

P23 

T 

P21' 

P21" 

P22’ 

P22” 

323’ 

P23” 

0.02 

0 

-0.04801 

-0.00094 

0 

0 

-0.000002 

.04 

0 

-.09602 

-.00376 

0 

0 

-.00002 

.06 

0 

-.l44o4 

-.00847 

0 

0 

-.00005 

.08 

0 

-.19208 

-.01506 

0 

0 

-.00012 

.10 

•  0 

-.24015 

-.02353 

0 

0 

-.00024 

.20 

0 

-.48101 

-.09413 

0 

0 

-.00191 

.30 

0 

-.72331 

-.21179 

0 

0 

-.00643 

.40 

0 

-,96774 

-.37652 

0 

0 

-.01525 

•  50 

0 

-1.21504 

-.58831 

0 

0 

-.02978 

.60 

0 

-1.46591 

-.84716 

0 

0 

-.05146 

.70 

0 

-1.72106 

-1.15308 

0 

0 

-.08171 

TABLE  II.-  VALUES  OF  KERNEL  AND  KERNEL  D 


MCA  TN  2213 


23 


d 

•\ 

oj  m  oovo  o  VD  fn-=J-  m  o  t—  o 
lAC^OOOQ^-t^OHHO  OnH 
t>-C>»OVO-3-  LTNCO  l>-0\t>-irN-^  CU 
IfNHOJt— HHITMTNOOOJH  OO 

ooc>str\ajHoooooo  oo 

HOOOOOOOOOOOO 

-^VDOh-OJOOOJCJNCVJOf-tt^ 

H  H  LTN  OvVO  OJOJ-?-t'-OVOCVJlA 

OJ  LTN  VO  H  O  LTN  t-"  LTN  on  t-VD  ON  C— 

O  O  H  oniTNCJN-?-  onH  on-?-  onia 
OOOOOOHOJ  on-?-  ITNVO  CO 

ooooooooooooo 

O  1  1  t  1 

1 

o 

CQ 

H  Os  VO  ITS  OJ  VO  COVD  OJ  VO  OO  ^-VO 

q  on  c^  o  ^  o\-^  H  onvo  h  vo  on 

H  H  lf\  0\-4-  ITNVO  ir\  H  OJ  O 

-=^O^oncvN^-OJOOO-:tOJHOO 
t—  ^-vo  -4-OOHHOOOOO 

ooooooooooooo 

o 

OJOO.OJ^.ITNWHOCVJLTNCVJA-?- 
CO  OJVO  AVO-?-  OJVO  At-VO  O  tH 

O  OJCOOOHVOH^VOCTN  t-VO  OJ 
o  o  o  H  onvo  H  ON  o  Q  on  OJ  on 
OOOOOOHH  on^  A  C—  CJn 

ooooooooooooo 

P 

O 

CQ 

H-=fOJ-=J-OHOJOJ  t— VO  CO  CO  VO 
f02S  ^  Q^co  Er-  ON  OJ 

irNco-=f  H  ONVo-=f  i>-ir\vo  oj 

P  O  ON^VOVOCOVO  OJITNVO  LfNO 
OOOITNOJ  O  ont-ONOJ  Ir-H 
HHHOOOOOOOHHOJ 

O  OJ  A  VO  O  OJoO  AVD  OJ -?-  A  H 

CO  .OJ  OJ  A  ON  CO  c—  t--  ON  H  t--  O 

HVD-?-AHOJCO-?'CJNAOOJ  on 

VO  A  VO  AOJVO  I>-OJH-?- 
HCO-?  OJH  CTNCOCOCOCOCO  ONO 
OJHHHHOOOOOOOH 

H 

O  1  1  1  1  1  1  1 

1  1  1  1  1  1  1  1  1  1  t  1  t 

t>- 

d 

U) 

p  t-vo  OJ  o  -=i-  ON  o 

P  CCJlTNONOJOJOJONlTNlfN 

pN^-COVO  HCO  OCO-=^VD  ^--s^  H 
t^ONOVO  COH  t— cn_=h^  ITNOnt— 
t>-ONOJ  cniAt^ONOJ  i>-oj  OJ  on-^ 
opt— JHiHHHOJOJ  on-^  H 

o  On-?-  t-^-?-  onvovD  o  oo -=t  ono 
d  gp  .CO  00  CVJ  o  H  on’co 

O  onvo  -?  OJ  o  OJ  O  A  0^1-?-  H 
i;;;-onA-?--?-  ononco  !^~-oncr^o  cvi 
-?-oncvjojt--cvjooJ-?-VDt>-ocvj 
ONAOJHOOOOOOOHH 

o  1  1  1  1  1  . . . 

1  1 

1  1  1  1  1  i  1 

t- 

• 

O 

•N 

CO 

onononoo  O  I>-OncO  onvo -:t  h vo 
Ono  O^ir-onoj  H  rooo^-t^o  H 
t^ONOniTNlTNOlNO  OLTNO  O  OJ  O 
-=!■  CO  ^-VD  lAcO  OJHcgonirNlTNO 
on-:?  ITN  LTN-?  OJ  O  OOVO  ON  OJ  t-  H 
OOOOOOOOOOHHOJ 

pvo  Acooo  OJco  onoo  h  on-?-  -?- 
VO  o  t-ACVJVOAOnHHoOA-?- 
cjn  H  t—  on  H  t—  H  ONVO  CO  on  on  A 

A  A-?-  CVJ  C--VO  onoo  O  O  -?-  t---  H 
HCO-?-  CVJOOO  D— AA-?-  onOJ  OJ 
OJHHHHOOOOOOOO 

o  1  1  1  1  1  1  1 

. . 1  1  1  1  1  1'  1 

t-- 

d 

CO 

^OOr^^-HOHOJ-=^ 
ITNNO  OOOJ  f-lTNlAO  OnVOoO  ONITN 
OOONONLfNLTNCOVDOOC^OVOHO 
oo  OnVO  on  IfNCO  CO  OO  ON  OJ  -=J-  cn  i>- 

on-?- 

OOOOiHHHOJOJ  m-?-  t--  H 

O  1  1  1  1  1  1  1  1  1  r  1* 

'  1 

^  ^  IH  H  ON^  AVD  .ON  OJ  CVJ  t^-VO 
ON^  n^-?--?•  OJ  t—ojVD  cnAOJO 

O  AAOn-?-  onoco  onAACO-ON 
^  o  00  S  ^  CVJ  aPS 

-=tmcvjcvjt--onooHcvjcvjojcvj 

C7NACVJHOOOOOOOOO 

1  1  1  1  1  1 

p 

CO 

ifN  t---  o  OJ  VO  ir\vo  VO  VO  On  m  o 

'P  ^  H  tTN  H  O  H 

onvo  o  t— -?■  ovovo  t— H  1-i  t-vo 

OJ^irVWonONOnONHLTNVD  HVO 

OJ  OJ  on-?-  mvo  onhvdon-?--?-  h 
ooqooooHHHOJ  on-?- 

P  t--  O  (TnVD  O  H  CO  00  t^--?-  CVJ  0\ 

Vp  O  CO  .ON  ONVO  A  ctnoo  On  on  b-oO 
^  on  A.VD  O  O  A  O  ACO  IH-?-  A 
^  A  onvo  ooojHooJONH-?-on 
CJNO  t— 00  H  C7NC0  OND— OJVOVO  CVJ 
ononHoooHOJ  on-?-  .?*.?--?■ 

w 

O  1  1  1  1  1  1  1  1  1  1  1  i" 

1 

1  1  1  1  1 

o 

CQ 

cornt-OJONtr-t^-H  ojvo  ojoo  m 
l:H^700JocQONH  on-?-  o  t--?- 
-?HonHOJ-^onocooNHont— 

H  OJ  on-?-  LfNOO-?-  OJOOlAONt-H 
OOOOOOHOJ  on  ITN  CO  ON-?- 

H  on  o  C7N-?-  -?-  C—  CVJ  A  H  D—  O 
'P  iDP  CVJ  CTsOnHH-  t- 

OJOOOOt— CTNOOOnoOHOOJVD 
AO\0nCJNO  CVIHA  ctn-?-  t--  t—  o 

H  vg  CJN  onvo  t— VO  on  h  o  co  vo  h 
On^  on  o  CO  vo  A-?-  on  CVJ  o  o  CVJ 

w 

o  H  on 

on  OJ  H  H  1  1  1  1  1  1  1 

1  1  1  1 

p 

d 

•\ 

CQ 

w 

-?-HOLrNVOO-?--?-ONOJOn-?-irN 
WVO  O  H  onirNt—lTNonONlTNO  OJ 
-?-  OJ  onoo  ONOO  OJVO  on  H  OJ  OnVO 
OJHOJonh-QVOOOOOrHVOVO 
HOJOJr?0-?-ONirN  onoo  t—  H  OJ 
oooooooHOJOJon  irwo 

O  1  1  1  1  1  t  1  1  1 

A  on  A  t— -?-  COCMAOVDHOOA 

OJ  H  A  A  CVJ  (On  OnVD  t—CO  A  H  -?- 

t>-.?-onocviojonHCJNOonHo 
q  qco  (ONAAco  hhctni^-d^-oj 
HONHOCVJOOonCVJoOOJono 
vo-?-oncvjHOHOjon  on-?-  ^  ^ 

1  1  1  1  1 

y— V 

E>- 

d 

•\ 

CQ 

P  ^  onv^  pjH-?-Lr\HNO  OVD  H 

VO  OOJQVO^-lfNO-?■COOJCOON 
ITNOJVO-?  ononOJt— OCO-?-  r?VO 

H  HVD  CVJ  t>-H-?-NO  OJCO-?  t^-cO 

S  d  irv  on  o  h  cm  onvo 

OOOOOOOOi~IOJ-?-OJCVJ 

A-?-vo  ONO  (jNCJNonoco  ono  m 
t-— OCO  O  ovovo  AVO  OJVO  OJVO 

H  A  OJ  t-  onvo  t-  H  CVJ  tH  A-?-  tH 
.co  A  t:-  CVJ  on  OJ  on-?-  vo  vo  oj  on  H 
vo  onvo  Hco-?-A-?  onoJ  h  mco 
CTNCJnH  ONt— vo  A-?-  onOJ  HO  H 

W 

Q  1  1  1  1  1  H  OJ 

CVJ  H  H  1  1  1  1  1  1  1  1 

1  1  I 

CO 

o  hc^h  ooNtr-oncT\H-?-  t--?- 

O  [>--?-  O  W  ITN-?-  CVJ  t--  H  irv  H 

pirNiHirNO-^ooH-?-irNt'~co  t- 
pqOt—OOOi-IOJt— HOOLTNOJon 
-?--^  lAonoJOco  oncjNH-?-  t--?- 
O  CVI  O  CO  ITN-?-  OJ  OJ  H  O  O 

^  HCJNont— CJNOHonHO 

o-w-?  A  w  o  t--?  P-o 
HCO-^  OAHAO 
on  pj  ago  Ht-OJHOOOtHCOQ 
d  o  CVJ  on  A^  -?^ 

OOHOJC\J-?-ACOOOJ-?-t--0 

OJi— IHHHl  1  1  1  1  1  1  1 

1  1  1  1  1 

H  H  H  H  CVI 

TABLE  IV.-  HINGE-MOMEET  COEFFICIENTS  FOR  FLAPPING  AND  ROTATION  OF  AIRFOIL 
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TABLE  V.-  VALUES  OF  COEFFICIEOTS  IN  SERIES  REPRESENTATION 

7.  =  -2V  ^0  2  ^  ^  ^  ^0^ 


[p„ = p„' + ip„”] 


For  flapping  of  airfoil 

(Dp  =  0.02 

(Op  =  0.04 

(Dp  =  0.06 

n 

n 

Pn" 

Pn' 

p  1  t 

P  * 

P  ^  * 

■^n 

0 

1 

2 

3 

k 

0.00405 

-.00050 

0 

0 

0 

0.02580 

.00005 

0 

0 

0 

0.01195 

-.00194 

0 

0 

0 

0.04726 

.00036 

.00003 

0 

0 

0.02102 

-.000413 

.00002 

0 

0 

0.06523 

.00091 

.00009 

0 

0 

CDp  =  0.08 

CD  =  0 . 10 

cDp  =  0.20 

(Dp  =  0.30 

n 

Pn' 

Pn" 

B 

Pn" 

Pn’ 

IBm 

Pn' 

■a 

0 

1 

2 

3 

4 

0.03053 

-.00700 

.00006 

0 

0 

0.08067 

.00181 

.00018 

0 

0 

0.03991 

-.01048 

.00013 

.00011 

0 

0.09420 

.00302 

.00034 

-.00004 

0 

0.08342 

-.03610 

.00123 

.00013 

0 

0.14570 

.01371 

.00218 

-.00007 

0 

0.12386 

-.07419 

.00454 

.00057 

-.00004 

0.18312 

.03356 

.00039 

-.00043 

-.00003 

05^  =  0.40 

COj.  =  0.50 

(Dj.  =  0.60 

(Dp  =  0.70 

n 

P  * 

P  ^  * 

Pn* 

p  f  t 

^  n 

Pn’ 

p  *  ^ 

P  ’ 

p  » * 

0.16431 

-.12273 

.01189 

.00167 

-.00018 

0.21075 

.06574 

.01353 

-.00157 

-.00014 

0.20492 

-.17824 

.02588 

.00378 

-.00059 

0.22786 

.11402 

.02355 

-.00432 

-.00044 

0.24175 

-.23186 

.04912 

.00634 

-.00165 

0.23334 

.17927 

.033JH 

-.01014 

-.00096 

0.27133 

-.27797 

.08434 

.00848 

-.00355 

0.22743 

.26203 

.03984 

-.02027 

-.00157 

MCA  TN  2213 


27 


TABLE  V.-  VALUES  OF  COEFFICIENTS  IN  SERIES  REPRESENTATION 


7  =  -2V^Pq  cot  2  ^  ^  Pn  -  Concluded 


For  rotation  of  airfoil 

CDp  =  0.02 

(Oj.  =  0.04 

=  0.06 

n 

p  ' 

P  *  ^ 

P  ' 

•^n 

p  ^  • 
rn 

Pn' 

P  t  T 

0 

1 

2 

3 

k 

1.29401 

.00272 

.00032 

0 

0 

-0.19220 

•03991 

-.00003 

0 

0 

1.19290 

.00793 

.00125 

0 

0 

-0.27984 

.07696 

-.00013 

0 

0 

1.10671 

.01391 

.00270 

0 

0 

-0.32806 

.11189 

-.00035 

0 

0 

(Dp  =  o.o8 

•  cOj.  =  0.10 

cDp  =  0.20 

coj.  =  0.30 

n 

p  * 

Pn" 

P  * 

-^n 

p  I  f 

P  ^ 

•^n 

P  '  * 

P  * 

■^n 

p  f  1 

■^n 

0 

1 

2 

3 

k 

1,03781 

.02012 

.00461 

0 

0 

-0.35794 

.14529 

-.00069 

-.00009 

0.97989 

.02691 

.00697 

-.00004 

0 

-0.37637 

.17790 

-.00126 

-.00018 

0 

0.80131 

.06175 

.02586 

-.00037 

-.00004 

-0.41701 

.33^05 

-.00612 

-.00135 

0 

0.70475 

.10655 

.05649 

-.00159 

-.00029 

-0.45044 

.48796 

-.01645 

-.00454 

,00014 

(%.  =  o.4o 

=  0.50 

CjDj,  =  0.60 

OD^  =  0.70 

n 

p  * 

P  '  * 

p  f 

■^n 

P  ^  * 

■^n 

P  * 

p  * ' 

P  * 

■^n 

P  *  ’ 

0 

1 

2 

3 

4 

0.63101 

.16912 

.09907 

-.00472 

-.00097 

-0.49399 

.64043 

-.03557 

-.01080 

.00055 

0.55644 

.25543 

.15194 

-.01190 

-.00196 

-0.54194 

.78449 

-.06873 

-.02105 

-.00169 

0,47233 

.36613 

.21279 

-.02383 

.00531 

-0.59052 

.91393 

-,11630 

-.03601 

.00378 

0.37484 

.49812 

.26823 

-.04688 

-.00929 

-0.61770 

1.00553 

-.19176 

-.05153 

.00848 
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Figure  1.-  Real  part  of  kgg  against  reduced  frequency  for  M  =  0.7. 
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Figure  11.-  Real  part  of  k. 


for  M  =  0.7. 
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Figur3  13.-  Real  part  of  k  for  M 

oR 


0.7. 
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